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Abstract 

The QCD corrections to electroweak parameters depend on the renormalization 
^ scheme and scales used to define the top-quark mass. We analyze these dependences 

for the IT-boson mass predicted via Ar to 0(aa s ) and 0{aa 2 s ) in the on-shell and MS 
schemes. These variations provide us with a hint on the magnitude of the unknown 
higher-order QCD effects and contribute to the theoretical error of the prediction. 
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-F; ' 1 Introduction 



• i 

It is well known that strong-interactions effects on the vacuum-polarization functions of the 
electroweak gauge bosons play a significant role in present and future high-precision tests of 
the standard model (SM) [l],§,[|(|]. In perturbative calculations to 0(aa s ) |5],||,|7||| , these 
effects arise from the type of two-loop diagrams where one virtual gluon is exchanged within 
a quark loop inserted into an electroweak-gauge-boson line. Since the top quark is by far the 
heaviest established elementary particle, with a pole mass of M t = (180 ± 12) GeV |§, the 
leading high-M 4 terms, of 0{G pM^ \ are particularly important. As for oblique corrections, 
i.e., those which arise from the gauge-boson vacuum polarizations, these terms together with 
their quantum-chromodynamical (QCD) corrections are all concentrated in Ap = 1 — 1/p, 
where p is the familiar parameter introduced in Ref. [|K|]. The leading-order QCD corrections 



to Ap have been known for several years MM; those of next-to-leading order have recently 
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been calculated ]TT| and found to be relatively large, having indeed a non-negligible impact 
on ongoing precision tests of the SM. It is of great phenomenological interest to estimate the 
residual theoretical uncertainty of the QCD corrections to Ap and other basic electroweak 



parameters such as Ar |fL2| . This is the motivation for the present paper. 

Q CD-improved analyses of electroweak parameters |],|^,|3],|^,|^,|],[7],[8] are usually carried 
out using the pole definition of the top-quark mass. This directly corresponds to the mass 
parameter which is presently being extracted with the Fermilab Tevatron f9j and will be 
with the CERN Large Hadron Collider (LHC) and future e + e~ linear colliders, since, in the 
propagation of the t and i quarks between their production and decay vertices, configurations 
near their mass shells are kinematically favoured. As a matter of principle, however, this 
mass convention is arbitrary, and we might as well adopt another one. For, if all orders of 
the perturbation expansion were taken into account, the final result should not depend on 
the selected scheme. Yet, this no longer holds true if the perturbation series is truncated. 
In general, the finite-order results also depend on the renormalization scales of the quark 
masses. It is generally believed that the scheme and typical scale variations may be used to 
estimate the theoretical uncertainty due to the unknown higher-order corrections. 

In this paper, we shall pay special attention to Ar which parameterizes the non- 
photonic corrections to the muon lifetime and allows us to indirectly determine the W- 
boson mass for given values of the top-quark and Higgs-boson masses. Using two popular 
definitions of quark mass in QCD, namely the pole mass and the mass of the modified 
minimal-subtraction (MS) scheme ||13|| , we shall quantitatively analyze the scheme and scale 
dependences of the Mw prediction at next-to-leading order in QCD, and so estimate the 
residual theoretical error on Mw from QCD sources. Our evaluation of Ar will make full 
use of the present knowledge of higher-order corrections, so that the extracted Mw values 
should be reliable within the quoted errors. Confrontation of these values with the future 
high-precision measurements of Mw at the Tevatron and the CERN Large Electron- Positron 
Collider (LEP2), in connection with a reduced error on M t , will allow us to pin down the 
mass of the SM Higgs boson and to facilitate the search for it. 



At this point, a few comments on the so-called tt threshold effects fT4j|Il] are i n order. The 
study of these effects was an attempt to estimate the dominant higher-order QCD corrections 
to Ap and other oblique electroweak parameters prior to their diagrammatical computation 
fTT| . This approach was based on the assumption that the bulk of the QCD corrections 
arises from the ladder diagrams of multi-gluon exchange in the tt system. The absorptive 
parts of these diagrams can be resummed in the non-relativistic approximation, producing 
a prominent enhancement of the tt excitation curve in the threshold region along with a 
lowering of its onset. This treatment naturally takes into account the finite lifetime of the 
top quark as well. The real parts of the diagrams were then found via dispersion relations 
with subtractions derived from Ward identities. The fact that the explicit 0{a 2 s GFMf) 



calculation of Ap ]TT[ nicely agrees with the rough tt-threshold estimate [T^] — in fact, it 
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comfortably lies within the errors quoted in Ref. |L4| — may be viewed as some posterior 



justification of this method. However, the good agreement is, to some extent, fortuitous, 
since approximately 30% of the O^oIGfM^) correction is due to double-triangle diagrams 
||1 1|| , which are beyond the scope of Ref. [|14]]. Therefore and for other reasons given in 
Ref. [|16|], we shall take the point of view that, at the present time, the result of Ref. JTT 



represents the most reliable description of the QCD corrections to Ap beyond the leading 
order, and that the residual theoretical uncertainty may be assessed by analyzing scheme 
and scale variations. 

This paper is organized as follows. In Section 2, we shall translate existing on-shell (OS) 
results for the gauge-boson vacuum polarizations in 0(aa s ) to the MS scheme of quark- 
mass renormalization. Furthermore, we shall take a look inside the renormalization-group 
(RG) structure of the QCD expansion of Ap to (9(a 2 GVM 2 ) [|TTJ. In Section 3, we shall 



quantitatively analyze the scheme and scale dependences of the M\y value predicted from the 
analysis of Ar to leading and next-to-leading order in QCD. Our conclusions concerning the 
theoretical uncertainty in My/ of QCD origin are summarized in Section 4. The Appendix 
contains some general relations which may be used to implement the scale dependence of the 
QCD coupling and mass in the MS scheme, and to switch between the MS and OS schemes. 



2 Formalism 

We shall work in the electroweak OS renormalization scheme, which uses the fine-structure 
constant, a, and the physical particle masses as basic parameters, and define (? w = 1 — = 
M^/Mf Pinl. Our analysis of Ar will be based on Refs. Following Refs. we 

shall use the decomposition 

c 2 

Ar = Aa - -f Ap(l - Aa) + Ar rem , (1) 

a w 

where Aa embodies the contributions from the charged leptons and the first five quark 
flavours that drive the fine-structure constant from the Thomson limit to the Z-boson scale, 
Ap is the top-quark-induced shift in the p parameter |TD[ written with the Fermi constant, 
Gf-i and Ar rem is devoid of large logarithmic and power-like terms of fermionic origin. We 
would like to point out that, unlike Ap, Ar rem must be written in terms of a, since it 
is via Ar that Gf is introduced into the SM |I|j2]j3|,p^1 . To our knowledge, all existing 
analyses to 0(aa s ) of Ar and other oblique electroweak parameters with a non-trivial quark- 
mass dependence, i.e., other than Ap, employ the pole definition of quark mass in QCD 
|l|J|J^,|,|,||,^,|| . In the following, we shall describe how these calculations may be converted 
to other schemes of mass renormalization in QCD, in particular to the MS scheme. 

The relevant quantities in this context are the transverse gauge-boson vacuum polariza- 
tions induced by a pair of quarks, with pole masses Mi and M 2 . Adopting the conventions 
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of Ref. 0, we may write the vector and axial- vector parts as 



n 



V,A, 



N, 



s ,M 1 ,M 2 ) = ^n 



NrC 



s, M h M 2 ) + a-^-?-nX' A (s, Mi, M 2 ) + 0(a 2 ), 



(2) 



where s is the square of the external four-momentum, N c = 3, Cp = (iV 2 — 1) /(2N C ) = 
4/3, and a = a s /7c. We shall employ dimensional regularization in n = 4 — 2e space-time 
dimensions and introduce a 't Hooft mass, /i, to keep the coupling constants dimensionless. 
We shall suppress terms containing 7^ — ln(47r), where 7^ is Euler's constant. These terms 
may be retrieved by substituting /i 2 — > 4ne~ lE /i 2 . In the MS scheme JIBfl , they are subtracted 
along with the poles in e. For later use, we list the lowest-order functions appearing in Eq. (|2|) 



47r 2 n^(s,M 1 ,M 2 ) 



e MiM 2 



+ 



1 + 



M 2 + Mf 



± 3MiM 2 



(Mf - M 
2^ 



2^ 2 



^ M 2 + M 2 -s 

V A arcosh 

2MiM 2 



Mf - M 



2s 

(Mf - M 



(Mf - M| 



I) 2 



2s 



2s 



+ 0(e) 



T 3M X M 2 



M 2 ^ 

In + -s - Mf - M 2 ± 6MxM 2 
M| 3 12 



(3) 



M 2 ^ 2 



where A = [s - (M x + M 2 f] [s - (Mi - M 2 ) 2 ]. Equation (§) is valid for s < (Mi - iw 2; 
It may be analytically continued to other values of s by observing that s comes with an 
infinitesimal, positive imaginary part, i.e., s — ► s + ie. Specifically, 



"\/A arcosh — - 



M| - s 



2M 1 M 2 



-A arccos 



Mf + Mf — s 
2MxM 2 



/- / s — Mf — Mf . \ 

-V A arcosh - - m \ . 

\ 2M X M 2 J 



(4) 



The general formulae for IlX' A (s, M 1: M 2 ) are somewhat lengthy ||. Fortunately, the finite- 
mass effects in 0(aa s ) on the oblique electroweak parameters (with the well-known exception 
of Aa) due to the first five quark flavours are exceedingly small []7|J§, so that, beyond 0(a), 
we may safely neglect all quark masses, except for M t . In this approximation, the mass- 
dependent cases may be cast into the form @ 



7T 



j^^Xis, M, M) = RX, + V l (R) + 0(e) 



7T 



-^n A ( S , M, M) = i?X x + n + Ai(it!) + (9(e) 



^i' A (s, M, 0) = \{XX l + Y l ) + F 1 (X) + 0(e) 



(5) 



where Vx, A x , and F x are finite functions of R = (s/AM 2 ) and X = s/M 2 , and |||T9 

X X 



1 + L _4C(3) + -, 
2e SV ; 12' 



2e 2 



3L 
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T 



n 
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3L 2 + — L + 6C(3) + -C(2) 
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(6) 



with L = ln(fi 2 /M 2 ). Here, £ is Riemann's zeta function, with values £(2) = 7r 2 /6 and 
C(3) « 1.202 057. 

The quantities in Eq. (|5]) refer to the OS scheme, i.e., they contain the contributions which 
emerge from the respective one-loop seed diagrams by inserting the OS mass counterterm, 
5M = m — M, where m is the bare quark mass, in all possible ways. The resulting shifts 



in U^ ,A may conveniently be constructed from Hq' A through the variation 



a ^cv jn p (s?Mi?M2 



N 2 d 

i=i l 



V,A, 



s,M 1 ,M 2 ). 



(7) 



For later use, we list the OS mass counterterm |20 
5M 



M 



—aCp 
—qCf 



3 3 



r 1 + e 



3-2e 
4e(l - 2e) 



+ 0(a 2 



46 



-L + 1 + e -L 2 + L 



-X(2) + 2)+0(e 



+ 0(a 2 



where T is Euler's Gamma function and L is defined below Eq. (||). 

In other mass renormalization schemes, m is differently split into the renormalized mass, 
m, and the counterterm, 5m, i.e., 



m = M + 5M = m + 5m, 



(9) 



where 5m differs from 5M by some finite amount. In general, m depends on the renormal- 
ization scale, /i, and on the QCD gauge parameter. Of all possible schemes, the MS scheme 
is singled out because there 5m just collects the e poles in Eq. (j|). Furthermore, the MS 
mass is gauge independent, as is M pij. We may translate Eq. to some other scheme by 
replacing Mj with m^ and adjusting the mass counterterms. The latter generates the shift 



N, 1 
3 



<, — - — AIL/ (s,m 1 ,m 2 ) = — 2_^{Mi - m;)-^— II ' (s,m 1 ,m 2 ), 



i=l 



(10) 
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where we have used Eq. (|9|). Notice that Eq. ( |T0| ) does not require knowledge of the 0(e) 



V,A 



terms of I\ 

In the following, we shall take m to be the MS mass. From Eq. 

'3. 



M = m 



1 + aC F [ -I + 1 



0(a 2 



we may read off that 
(11) 



where / = ln(/i 2 /m 2 ). The 0(a 2 ) term of Eq. (pTTj) is written out in Eq. fl3"0|). With the help 
of Eqs. (H), (0), and (0), we may now convert Eq. (|) to the MS scheme. The resulting 
expressions emerge from Eq. (|5|) by replacing M, R, X, X\, Y%, Vi(R), Ai(R), and Fi(X) 
with m, r = (s/4m 2 ), x = s/m 2 , 



1 

X 1 = - + i-4 C (3) + -, 



75 



^i(r) = ^(r) + (3/ + 4) 



arcsin A /r 



+ 1 



rJl/r — 1 



A 1 (r) = A 1 (r) + (3Z + 4) 24 



Fi(x) = Fi(z) + -(3Z + 4) 



x 



1 arcsin i/r — 1 

1\, , , I 
— J ln(l -x) 



x 



12) 



respectively, where / is defined below Eq. (|TT]). The last three lines of Eq. (|T^) are valid for 
< r, x < 1, the range relevant for the evaluation of electroweak parameters. Expressions 
appropriate for other values of r and x may be obtained with ease by analytic continuation 
as described in Ref. The following special cases frequently occur in applications: 



Vi(0) = V1(0) =0, 
7i(0) = ^(0)-|(3/ + 4), 
A 1 (0)=A 1 (0) + 3l + 4, 
F 1 (0) = F 1 (0) + i(3Z + 4). 



(13) 



In the remainder of this section, we shall discuss the incorporation of three-loop QCD 
corrections |TT|j22]| in Eq. ([!]). In Ref. [|ll|], the 0(a 2 GpM 2 ) term of Ap has been calculated 
using both the OS and MS definitions of quark mass in QCD. In Ref. [p2]] , the 0(aa 2 ) 
correction to Ar> em has been expanded in powers of M|/M 2 , and the first two terms of this 
expansion have been presented both in the OS and MS schemes. It is instructive to rewrite 
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the results of Ref. [|TTJ in such a way that the RG origin of the logarithmic terms as well 
as the interrelation between the OS and MS versions are explicitly displayed. This may be 
achieved with the help of Eqs. (|22|) and (^) for rif — 6 active quark flavours and leads to 



Ap = N c X t l + aR 1 (l + a(3 L) + a 2 R 2 + 0(a 3 



with 



3X t 



I - 2.8.->9 9.12r/ ( 1 + -aL 



14.594 028 a 2 + 0(a 3 



3x t H + a 



21 - 0.193 245 1 
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-al 



(14) 



Ap = N c x t {l + a [2 7o Z + n (l + a(p + 2 7o )/) 
+ 0(a 3 )} 



7o (A) + 2 7o )Z 2 + 2(-2 7o 2 + 7l )Z + r 2 



15 
T 



Z 2 + ^Z- 3.969 560) + 0(a 



(15) 



r x = i?i + 2C F , 

r 2 = R 2 + C F {2R 1 - 4 7o + C F ) + 2iT , 



(16) 



where X t = (G F M? /8n 2 V2), x t = (G F m 2 t /8tt 2 V2), and L and I are defined below Eqs. (§) 
and (pH), respectively. The first three coefficients of the beta function and the quark- mass 
anomalous dimension of QCD are listed in Eqs. ( p0|) and (pl|), respectively, and is defined 
in Eq. (p9|) . The genuine information on the QCD correction to Ap is carried by R\ [||,[7| 
and i? 2 [|TT|] or, equivalently, by r\ [^] and r 2 [|TT| . 

In the OS scheme, all p dependence originates in the renormalization of a, while the MS 
scheme has m as an additional source of p dependence. So far, we have used a common 
renormalization scale, p, for a and m. In the following, we shall abandon this restriction 
and distinguish between coupling and mass renormalization scales, p c and p m . We shall 
henceforth use the symbol p c in the OS expressions. In order to disentangle p c and p m in 
the MS formulae, we first replace p with p m and then substitute 



a(p m ) = a(p c ) 



l-a(p c )/3 Q ln% + 0(a 2 



(17) 



which follows from Eq. fl2"2"|). Finally, we expand the resulting expressions in powers of a(p c ) 
and truncate these expansions beyond the order under consideration. 

In the next section, we shall estimate the QCD-related theoretical uncertainty in My/ by 
studying the scheme and scale dependences of the Ar analysis. In order for our estimate 
to be meaningful, we shall have to judiciously choose the central values and widths of the 
scale intervals. It is natural to define the central values of p c and p m in such a way that 
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there the radiative corrections are devoid of logarithmic terms. Looking at Eqs. fll4"|) and 
(fjjD, we are thus led to set /i c = £, c M t in the OS scheme, and /i c = ^ c fi t and /i m = ^mfit, 
where fx t = mt{U-t), in the MS scheme. Here, £ c and £ m are variable numbers of order unity. 
A closed expression for fl t in terms of M t may be found in Eq. fl31~|). 



3 Numerical analysis 

We are now in a position to quantitatively explore the scheme and scale dependences of 
the M w prediction on the basis of Ar. Our starting point is the OS analysis of Ref. 0, 
which includes the dominant corrections beyond one loop, of 0(aa s ) f7j, 0{G 2 F Mf) [p4j] , 



and O^GpM^Mfj) fI],P|p5|| . Here, we extend this analysis to 0(aa 2 ) by accommodating the 
respective corrections to Ap [|Il|] and Ar rem |22|] in Eq. ([J). Furthermore, we update our 



input parameters according to Refs. PB"| , |2"T| , P8| , P9"] . In particular, we use the combined LEP1 



value M z = 91.1887 GeV |27|, the value Aa^ d = 0.0280 for the hadronic contribution to 



Aa, and the world average af\M z ) = 0.118 |p9fl . For the purpose of studying the scheme 
and scale dependences of Ar, we may assume that M#, M t , and a^'(Mt) are precisely 



known. Unless stated otherwise, we choose Mh = 300 GeV |27| and M t = 180 GeV ||. 
We compute af\M t ) at three loops in two steps. First, we scale af\iJ, c ) from p, c = Mz to 
/i c = M t via Eq. (|23|) with rif = 5. Then, we cross the flavour threshold at fi c = M t using 
the matching condition fl32j). For M t = 180 GeV, we so obtain af\M t ) = 0.1071. In order 
for our findings concerning the scheme and scale dependences to be meaningful, it is crucial 
that we consistently evaluate a^(n), m t (fx), and fx t to the order that we consider at a time. 
For the analysis of Ar to 0(aal), we first compute a^(/i) using Eq. ( p3[ ) forn/ = 6 with the 
O(a^) terms in the denominator omitted. Then, we insert this value together with a^(M t ) 
and m t (M t ), which we extract from Eq. (]2"8|), into Eq. fl2"5D to obtain m t (n). We evaluate jit 
by means of Eq. ([H]). If we calculate Ar to 0(aa s ), then we correspondingly discard the 
0(a 2 s ) terms in Eqs. ([23]), and (0) and employ Eq. (0) instead of Eq. 

In Table p], we investigate the scheme dependence of the M\y prediction via Ar to 0(aa s ) 
and 0(aa 2 s ) within the ranges 160 GeV < M t < 210 GeV and 60 GeV < M H < 1 TeV. For 
the time being, we suppress the logarithmic terms of RG origin by choosing £ c = £ m = 1. For 
this choice, the MS values are always in excess of the the respective OS numbers, i.e., the 
QCD corrections are less negative in the MS scheme. In Fig. [I], we display this excess as a 
function of M H for M t = 168, 180, and 192 GeV g both in 0(aa s ) and 0(aa 2 s ). We observe 
that the scheme dependence increases both with M t and Mh- For the M t and Mh intervals 
considered in Table [l], it ranges from 4.7 to 13.2 MeV in 0(aa s ) and from 2.2 to 10.3 MeV 
in 0(aa 2 ). As expected, it is significantly reduced as we pass from 0(aa s ) to 0(aa 2 ). 
For M t = 180 GeV, the reduction amounts to 51% (28%) for M H = 60 GeV (1 TeV). In 
Table [[], we also list the values of Ji t m 0(a s ) and 0(a 2 ). In 0(a 2 ), they lie by a roughly 
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Table 1: Mw predicted for various values of M t and Mh via Ar to 0(aa s ) and 0(aa 2 ) 
in the OS and MS schemes with £ c = £ TO = 1. For reference, also the Mw values evaluated 
without perturbative QCD corrections are given. Also p, t is listed. All masses are given in 
GeV. 





Mt 


M H 


M w 


0(a s ) 


0(a 2 ) 


w/o QCD 


0(a s ) OS 


0(a s ) MS 


0(a 2 s ) OS 


0(a 2 s ) MS 


160 


153.1 


151.5 


60 


80.398 


80.343 


80.348 


80.332 


80.335 








300 


80.295 


80.240 


80.245 


80.229 


80.232 








1000 


80.199 


80.143 


80.150 


80.133 


80.137 


170 


162.7 


161.1 


60 


80.468 


80.407 


80.413 


80.396 


80.398 








300 


on o£?o 

80.362 


80.302 


on onn 

80.309 


on onn 

80.290 


80.294 








1000 


80.265 


80.204 


80.212 


80.193 


80.198 


180 


172.3 


170.6 


60 


80.540 


80.474 


80.480 


80.462 


80.464 








300 


80.433 


80.367 


80.374 


80.354 


80.359 








1000 


80.334 


80.268 


80.276 


80.255 


80.261 


190 


182.0 


180.2 


60 


80.616 


80.544 


80.550 


80.530 


80.534 








300 


80.506 


80.434 


80.443 


80.420 


80.426 








1000 


80.405 


80.333 


80.343 


80.319 


80.327 


200 


191.6 


189.7 


60 


80.695 


80.616 


80.624 


80.602 


80.605 








300 


80.582 


80.504 


80.513 


80.489 


80.496 








1000 


80.479 


80.400 


80.412 


80.386 


80.395 


210 


201.2 


199.3 


60 


80.777 


80.692 


80.700 


80.676 


80.680 








300 


80.661 


80.576 


80.587 


80.560 


80.568 








1000 


80.555 


80.470 


80.483 


80.454 


80.465 



constant amount of 10 GeV below the respective M t values. This may also be seen from 
Eq. (|3l"l) , where, in the case of top, the coefficients of a s (M t )/ir and [a s (M t )/ii] 2 are —4/3 
and —6.458 784, respectively. To assess the significance of the scheme and scale dependences 
relative to the overall QCD effect on Mw, we also include in Table [l] the respective Mw 
values evaluated with the (perturbative) QCD corrections switched off. Of course, we cannot 

(5) 

reliably eliminate the intrinsic non-perturbative QCD corrections contained in Aa h( ^ d , i.e., 
without introducing a significant dependence on the poorly-known light-quark masses. 

In the remainder of this section, we shall stick to the central values M t = 180 GeV and 
M H = 300 GeV, and consider scale variations with 1/8 < £ c ,£ m < 8. In Table |2] and Fig. ^ 
we study the [i c dependence of Mw to 0(aa s ) and 0(aa 2 ) in the OS scheme. Obviously, 
the scale dependence is dramatically reduced, from 33.7 MeV to 11.1 MeV, i.e. by 67%, if 
the 0(aa 2 ) correction is taken into account. This should be compared with the shifts in Mw 
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Table 2: M w (in GeV) predicted for M t = 180 GeV and M H = 300 GeV via Ar to 0(aa s ) 
and 0(aa 2 ) in the OS scheme with £ c variable. 







GeV] 


e>K) 


0{al) 


1/8 


80.346 


80.351 


1/4 


80.354 


80.351 


1/2 


80.361 


80.352 


1 


80.367 


80.354 


2 


80.372 


80.357 


4 


80.376 


80.359 


8 


80.379 


80.362 



induced by the QCD corrections to 0(aa s ) and 0(aa 2 ) for \i c = M t , which are —66.1 MeV 
and —78.6 MeV, respectively, as may be seen from Table [l|. It is interesting to note that 
2.4 MeV, i.e., 19%, of the difference between the 0(aa s ) and (D(aa 2 ) evaluations of Mw 
is due to the three-loop correction to Ar rem The 0(aa 2 ) evaluation exhibits a local 

minimum at £ c = 0.176. This is the point advocated by the principle of minimal sensitivity 
(PMS) pD| . The 0(aa s ) and 0(aa 2 ) curves cross over at £ c = 0.183, the point of fastest 
apparent convergence (FAC) |31| . In the OS analysis of Ap to 0(a 2 GpM 2 ) with rif = 5, 
these points occur at £ c = 0.224 and £ c = 0.264, respectively P2"| . We note in passing that 
the application of the Brodsky-Lepage-Mackenzie (BLM) scale-setting criterion to Ap 
leads to £ c = 0.154 []lT| , |32| , p4| , which had been anticipated in the pioneering work of Ref. p4| 
prior to the advent of the 0(a 2 GpM 2 ) calculation of Ap | ]TT| . 

In Table |3| and Fig. ||, we investigate how the 0(aa s ) and 0(aa 2 ) calculations of M w 
in the MS scheme depend on p c and p m . We notice that the p c dependence is rather feeble 
for fi m fa fl t . This may be understood by observing that the coefficient of a s (fi c )/ir in the 
QCD expansion of Ap in Eq. (|1^) is then greatly suppressed f2~3fl . In Fig. |3], the points of 
minimal sensitivity, i.e., with zero tangents, are marked with "x." They are saddle points 
and gathered in a small strip around p m = Jx t . Their (£ c , £ m ) coordinates are (2.732, 1.024) in 
0(aa s ) and (0.133,0.853) and (7.127, 1.147) in 0(aa 2 ). For fixed £ c < 1, the 0{aa s ) value 
of M w varies quite strongly with £ m , by 185.4 MeV for £ c = 1/8. This has to be compared 
with the shift in Mw due to the 0(aa s ) correction for £ c = £ m = 1, which only is 58.6 MeV 
in size (see Table |l|). Of course, such an extreme variation cannot be interpreted as the 
uncertainty due to the neglect of higher-order QCD corrections. This rather tells us that 
our choice of the £ m interval width is not judicious in this case. Fortunately, the 0(aa 2 ) 
calculation is much more stable under scale variations. Here, the overall fluctuation is just 
25.7 MeV, while the QCD-induced shift in M w for £ c = £ m = 1 is 73.8 MeV in magnitude 
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Table 3: M w (in GeV) predicted for M t = 180 GeV and M H = 300 GeV via Ar to 0(aa s ) 
(upper entries) and 0(aa 2 ) (lower entries) in the MS scheme with £ c and £ m variable. 





1/8 


1/4 


1/2 


1 


2 


4 


8 


1/8 

I/O 


80 256 


80 304 

UU . kj Utc 


80 343 


80 374 

UU . ' ) I t: 


80 399 

UU >KJ u u 


80.421 


80 441 

UU . tit: _L 




80.365 


80.360 


80.357 


80.357 


80.357 


80.357 


80.357 


1/4 


80 200 


80 325 


80 352 


80 374 

UU .<J 1 t: 


80 392 


80 407 

UU . t:U 1 


80.421 




80.349 


80.352 


80.355 


80.358 


80.361 


80.363 


80.366 


1/2 


80.318 


80.342 


80.360 


80.374 


80.386 


80.395 


80.404 




80.343 


80.349 


80.354 


80.358 


80.362 


80.366 


80.369 


1 


80.342 


80.356 


80.366 


80.374 


80.380 


80.386 


80.390 




80.344 


80.350 


80.355 


80.359 


80.363 


80.366 


80.368 


2 


80.362 


80.368 


80.372 


80.374 


80.376 


80.377 


80.378 




80.347 


80.352 


80.356 


80.359 


80.362 


80.364 


80.366 


4 


80.379 


80.378 


80.377 


80.375 


80.372 


80.370 


80.369 




80.353 


80.356 


80.358 


80.360 


80.361 


80.362 


80.363 


8 


80.394 


80.387 


80.381 


80.375 


80.369 


80.364 


80.360 




80.360 


80.360 


80.360 


80.360 


80.360 


80.360 


80.359 



(see Table p. However, this variation is still approximately 2.3 times as large as the one in 
the corresponding OS calculation. 

Since it is hard to extract precise numbers from the contour plots in Fig. [|, we list in 
Table [| the maximum deviations of the MS evaluations of Mw to 0(aa s ) and 0(aa 2 ) within 
the variable range l/£ ma x < £o£m < imax from the respective values for £ c = £ m = 1. For 
completeness, we also list the corresponding numbers for the OS analyses to 0(aa s ) and 
0(aa 2 s ) of Fig. [| It is interesting to observe that, before the appearance of the 0(aa 2 ) 
corrections in Eq. (H) [[y])il]> the uncertainty due to the lack of these terms could have been 



estimated from the scale variation of the 0(aa s ) OS (MS) calculation with £ maa; = 4 (2.1). 
Thus, we may expect that similar scale variations will also yield meaningful results in the 
next order. 

The separation of /x c and \x m is perhaps not so easy to motivate on physical grounds. In 
Fig. |], we analyze the scale dependence of the MS calculation of M w to 0(aa s ) and 0(aa 2 ) 
identifying \i c = \i m = ^Ji t . We observe that the 0(aa s ) analysis is very unstable for £ <C 1. 
For £ = 1/8, we have M w = 80.256 GeV (see Table |3|), which is way below the M w range 
considered in Fig. [|. On the other hand, the 0(aa 2 ) curve nicely oscillates around the M w 
value at £ = 1, with a band-width of 13.5 MeV. In this one-dimensional analysis, the PMS 
points appear at £ = 1.601 in (D(aa s ) and at £ = 0.280 and £ = 2.898 in 0(aa 2 ). They are 
indicated by "o" in the contour plots of Fig. [| There is one point of FAC within Fig. |], at 
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Table 4: Maximum deviations (in MeV) of the My/ values predicted for M t = 180 GeV and 
Mh = 300 GeV via Ar to 0(aa s ) and 0{aa 2 s ) in the OS and MS schemes with l/^ max < 
£« < imax from the respective values for £ c = £ m = 1. 



lo§2 £,max 




AM W 


[MeV] 






C(a s ) OS 


<D(a s ) MS 


0{a 2 s ) OS 


0{otf) MS 


1/2 


+2.5 


+4.5 


+1.2 


+1.9 




-2.8 


-5.2 


-1.1 


-2.4 


1 


+4.9 


+ 11.4 


+2.4 


+3.6 




-5.8 


-14.4 


-2.1 


-5.0 


3/2 


+7.0 


+20.8 


+3.7 


+5.3 




-9.0 


-28.6 


-2.9 


-7.4 


2 


+9.1 


+32.9 


+5.0 


+6.9 




-12.6 


-49.4 


-3.4 


-10.0 


5/2 


+ 10.9 


+48.0 


+6.2 


+8.4 




-16.6 


-78.6 


-3.6 


-12.8 


3 


+ 12.7 


+66.5 


+7.5 


+9.8 




-21.0 


-118.8 


-3.6 


-16.0 



£ = 0.413. Another one is located at £ = 8.358. The BLM criterion only applies to /i c , but 
not to fi m , so that its implementation in our MS analysis is ambiguous |35] . 



4 Conclusions 

In this paper, we have extended an existing calculation of Ar in the OS scheme to 

We have then 



next-to-leading order in QCD by incorporating new three-loop results |Tl|j2 
converted this analysis to the MS scheme of quark-mass renormalization in QCD. Armed 
with these results, we have analyzed the scheme and scale dependences of the M w value 
predicted for given values of M t (pole mass) and Mh- We have verified that both scheme 
and scale dependences are considerably reduced if the next-to-leading-order QCD corrections 
are taken into account. For M t = 180 GeV and Mh = 300 GeV, the scheme dependence 
at the central renormalization point, £ c = £ m = 1, is reduced from 7.5 MeV to 4.8 MeV 
(see Table |T] and Fig. [T]). The scale dependences within the interval 1/4 < £ c ,£ m < 4 are 
decreased from 21.7 MeV to 8.4 MeV in the OS scheme and from 82.2 MeV to 16.9 MeV in 
the MS scheme (see Tables |, § | and Figs, g |). 

Although the numbers presented here are uniquely determined by the generally accepted 
rules of perturbation theory, their interpretation in terms of a QCD-related theoretical un- 
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certainty in Mw allows for a certain amount freedom. Nevertheless, we shall propose an 
algorithm to extract a central value and an error of Mw- If we adopt the point of view that 
the OS and MS analyses represent two independent theoretical determinations of M w with 
their individual errors, then we may combine them assuming Gaussian statistics as we usu- 
ally do with independent experimental measurements |26||. It is plausible and conservative 



to identify the error in each scheme with the absolute of the maximum deviation of Mw in 
the interval l/£, ma x < £ c , £m — imax from the central value, at £ c = £ m = 1. This will leave 
room for the reader to select his preferred value of £ max . For illustration, we shall again 
assume that M t = 180 GeV and M H = 300 GeV. The OS and MS central values of M w 
and the respective errors may then be read off from Tables |2], |3], and f|. For £ max = 4, we so 
obtain (80.367 ± 0.012) GeV to 0(aa s ) and (80.355 ± 0.004) GeV to 0(aa 2 s ). We would like 
to point out that the 0(aa 2 s ) central value is encompassed by the 0(aa s ) error, and that, 
in each order, the OS and MS central values lie within the error. This reassures us of the 
soundness of our procedure. The corresponding results for other values of £, max may be seen 
from Fig. |5|. It is clear that this approach does not make sense if £ max is so small that the 
scale variations of the OS and MS calculations do not overlap. 

Apart from the theoretical uncertainty due to the lack of knowledge of QCD corrections 
to Ar beyond three loops, which we have estimated here, there is another source of QCD- 
related error, namely the one connected with Aa^ dJ which is extracted from experimental 



data of e + e — * hadrons via a dispersion relation p8],[3(|. In Ref. [p8| , this error has been 
estimated to be ±0.0007, which, for M t = (180 ± 12) GeV and 60 GeV < M H < 1 TeV, 
translates into an error of approximately ±13 MeV in Mw- Other recent determinations of 
^ a hld [0 agree with the result of Ref. within less than two standard deviations of the 
latter. 

Finally, we would like to mention that the lack of knowledge of three-loop and subleading 
two-loop electroweak corrections to Ar represents another source of theoretical error on the 
Mw determination from the measured muon lifetime. The magnitude of the recently calcu- 
lated 0(G 2 F M 2 M w ) correction to Ap |57| indicates that the still uncontrolled 0(G 2 F M 2 M w ) 
term of Ar might have the potential to jeopardize the accuracy of the indirect Mw determi- 
nation. The study of the scheme and scale dependences of the dominant two-loop electroweak 



corrections provides us with a clue to the size of higher-order electroweak effects [38 
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A Appendix 

In this Appendix, we shall provide a few general relations which are useful for implementing 
the \x dependence of the QCD coupling constant and the quark masses in the MS scheme. 
We shall take the colour gauge group to be SU(iV c ); C F = (N 2 — 1)/(2N C ) and Ca = N c 
are the Casimir operators of its fundamental and adjoint representations, respectively, and 
Tp = 1/2 is the trace normalization of its fundamental representation. We shall keep the 
number of active quark flavours, rif, arbitrary. The RG equations for the so-called couplant, 
a(/x) = a s (/i)/7r, and m(/i) may be found, e.g., in Ref. [39|. For the reader's convenience, we 
list them here: 



where f|0] 

A> = 
A = 

02 = 



1 

4 
1 

16 
1 

64 

+ 



da 
din /i 
dlnm 
dlnfi 



[3(a) = -a 2 [p + foa + (3 2 a 2 + 0(a 3 
-7m(a) = -a 7o + 71a + 7 2 a 2 + C(a 3 ) 
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C A - \ T Fn f ^j 



34 



20 



AC F T F n f - —C A T F n f ) 



(18) 
(19) 



— C! 2 
3 Ca 

2 ^C\ + 2C 2 T F n f - ?fc F C A T F n f - l -^C 2 A T F n f + fc F T 2 n 2 



158 



C A T 2 F n) 



27 



(20) 



are the first three coefficients of the Callan-Symanzik beta function and |20| , f4T 



7o 
7i 

72 



97 



10 



16V2 F 6 3 1 

1 r 129 

64 l~2~ 



129 



C f Ca + 



11413 



108 



C F C\ + C 2 F T F n f {48((3) - 46] 



+C F CAT F nf 



-48C(3) - ^ 
sv 1 27 



140 2 2 \ 



(21) 
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are the first three coefficients of the quark-mass anomalous dimension. Here, ( is Riemann's 
zeta function, with value £(3) ~ 1.202 057. Our aim is to evaluate a = a(/x) and m = m(/x) 
for /x arbitrary, assuming that their values, ao = a(/xo) and itlq = m(/xo), are known at some 
starting scale, /x . 

Beyond the leading order, Eq. (|T8| ) cannot be exactly solved for a. A perturbative solution 
reads 



a = a 



1 - aM + a 2 £([3 2 £ - (3,) + a 3 £ (-/? 3 £ 2 + ~/W - &) + O(a 4 £ 4 ) 



(22) 



where i = ln(/x 2 / /Xq) . The leading logarithms in Eq. (p2|) may be resummed by writing [52 

a 



1 + a £ \/3o + a A + a§(-/W/2 + &)] ' 



(23) 



We recover the exact leading-order solution of Eq. (|1|) by discarding the terms of O{ofy in 
the denominator of Eq. (^) . 

Knowing the /x dependence of a, it is sufficient to obtain m as a function of a. Dividing 
Eq. ([19]) by Eq. fllBf ), we obtain a differential equation which may be exactly solved for m, 
given the coefficients of (3 and 7 m to a certain order. The exact solutions to leading and 
next-to-leading orders read 



m = m \ — 

M 



7o /A) 



m = m 



a 



tf + oqPx 



(24) 
(25) 



respectively. A perturbative solution of Eq. (TT9|) similar to Eq. (^) is given by 



m = m I 1 — a 7o^ + a 2 / 



7o 



(A) + 7oK - 7i 



+ a 3 l 



„ (PI , A,7o , 7 2 ' 



+ ( Am + ^7? + 7o7i ] ^ - 72 



+ W 4 ) 



(26) 



where £ is defined below Eq. ( |22"|) . By iterating Eq. ( p6|) , we may generate a closed expression 
for the mass parameter p, = m(/x) in terms of a and m for /x arbitrary: 



/x = m < 1 + 070/ + a / 



y(A) + 7o)/-27o 2 + 7i 
#l7o 



+ a 3 / 



+ -3/?o7 + A)7l + ^f" - 2 7o + 7o7i I + 47 - 47o7i + 72 



+ 0(a 4 / 4 ) , (27) 
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where I = ln(/i 2 /m 2 ). The right-hand side of Eq. (|27| ) is manifestly RG invariant through 
0(a 3 ). 

In the remainder of this section, we shall consider QCD with one massive quark and 
rif — 1 massless flavours. The relation between the pole mass, M, and m(M) at next-to- 
leading order is given by p]|.[l3[ 



M 



m(M) 



1 + a(M)C F + a 2 (M)K + 0{a 3 ), 



(28) 



where 43 



K = C F 



|C(2) - ? 
4 SV ; 8 



Cl 



|C(3) + C(2) ("31n2 + ^ 



121 



128 



|c(3)+«2)(|ln2-I 



+CfCa 
17.151430- 1.041 367 n f 



1111 
~3M 



C F Uf 



C(2) 71 



4 



192 



(29) 



Here, £(2) = tt 2 /6 and the value of C(3) is listed below Eq. (^). Using Eq. ( |28| ) along with 
Eqs. ([22] ) and (|26|) , we may express M entirely in terms of MS parameters, viz. 



M = m\l + a(j Q l + C F ) + a 2 



y Wo + lo)l 2 + (-27o + 7i + C F ([3 + 7o )) Z - 2 7o Cf 



3/3n 



(30) 



where I is defined below Eq. (^). The right-hand side of Eq. (|30|) is RG invariant through 
O (a 2 ) as it must be. Substituting Eq. (|28D into Eq. ( ^7|) evaluated at = M, we obtain a 
closed expression of ft in terms of M [fHjj , 

(31) 

In Eqs. (^8|), fl30|), and (|3i"D, a refers to rif quark flavours, i.e., a = a^ nf \ Finally, we list 
the matching condition for a at the heavy-flavour threshold, /i = M. It follows from the 
decoupling relation found in Ref. [^] and reads 



fl = M {l - a(M)C F + a 2 (M) [C F (2j + C F ) - K ] + C(a 3 )} . 



t W)(M) = a (n '~ 1) (M){l + 



^(M^T^^-^l+^a 3 ^ 



(32) 
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FIGURE CAPTIONS 



Figure 1: Difference (in MeV) between the MS and OS evaluations of Mw via Ar to 0(aa s ) 
(dashed lines) and 0(aa 2 ) (solid lines) with £ c = £ m = 1 as a function of M H for M t = 168, 
180, and 192 GeV. Upper curves correspond to larger M t values. 

Figure 2: /i c dependence of the OS evaluation of M w (in GeV) via Ar to 0(aa s ) (dashed 
line) and O(aaf) (solid line) for M t = 180 GeV and M H = 300 GeV. 

Figure 3: fi c and /i m dependences of the MS evaluation of Mw via Ar to (a) 0(aa s ) and 
(b) 0(aa 2 ) for M t = 180 GeV and Mh = 300 GeV. The contours of constant deviation 
(in MeV) from the value at £ c = £ m = 1 (marked by +) are shown in the (log 2 £ c , log 2 £ m ) 
plane. The saddle points in this plane are marked with the maxima and minima on the 
diagonal ^ c = £ m are marked with "o." 

Figure 4: /j, c dependence of the MS evaluation of M w (in GeV) via Ar to 0(aa s ) (dashed 
line) and 0{aa 2 s ) (solid line) with // c = fx m for M t = 180 GeV and M H = 300 GeV. 

Figure 5: Central values and QCD-related errors of M w (in GeV) evaluated via Ar to 
0(aa s ) (dashed lines) and 0(aa 2 s ) (solid lines) with l/£, max < £ c ,£ m < £ ma x for M t = 
180 GeV and M H = 300 GeV as a function of £ rnax . 
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